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Stochastic exclusion processes play an integral role in the physics of non-equilibrium statistical mechanics. 
These models are Markovian processes, described by a classical master equation. In this paper a quantum 
mechanical version of a stochastic hopping process in one dimension is formulated in terms of a quantum master 
equation. This allows the investigation of coherent and stochastic evolution in the same formal framework. 
The focus lies on the non-equilibrium steady state. Two stochastic model systems are considered, the totally 
asymmetric exclusion process and the fully symmetric exclusion process. The steady state transport properties 
of these models is compared to the case with additional coherent evolution, generated by the X X -Hamiltonian. 

PACS numbers: 03.65 Aa, 03.65.Yz, 05.50.Gg, 47.70.Nd 



I. INTRODUCTION 



Stochastic exclusion processes have been studied in statis- 
tical mechanics for a long time (HE). These are simplified 
one-dimensional hopping models, that allow the study of non- 
equilibrium phenomena in many -particle systems. The ex- 
clusion processes are modeled by a classical master equation 
that determines the time evolution of the probability distri- 
bution. The steady state of the master equation exhibits in- 
teresting non-equilibrium behavior, such as the presence of a 
current, non-equilibrium phase transitions and entire phases 
with a diverging correlation lengthen). The presence of cur- 
rents, such as the current of particles, energy or momentum, 
is a common feature of non-equilibrium steady states and can 
have profound effects on the correlations present in the sys- 
tem (6] [3 . Non-equilibrium systems can develop long-range 
correlations in the presence of a high current. 
The asymmetric exclusion process (ASEP) as well as the sym- 
metric exclusion process (SEP) are prime examples for such 
model systems EJEJ. These processes describe the hopping of 
hard-core particles in a one-dimensional chain, only driven by 
the inflow and the outflow of particles at the boundaries of the 
chain. Here one considers open boundary conditions, where 
particles are injected at the first site and are removed at the 
last site N of the chain. The steady state properties are en- 
tirely determined by the inflow and outflow. The dynamics of 
the particles in the bulk are given by translationally invariant 
hopping rates, that either constrain the hopping of particles to 
take place in only one direction (ASEP) or allow for a hopping 
in both directions (SEP). 

Transport properties of open quantum mechanical systems, on 
the other hand, are subject to recent research activities. A 
general interest is placed on how external noise, generated by 
the environment, affects the coherent transport in the system. 
It has been found, that the presence of noise in the quantum 
mechanical systems can actually aid the transport process of 
excitations through heterogeneous environments GOO, such 



as bio-molecules. An optimal ratio between coherent trans- 
port and dephasing noise can be found. The dynamics of 
open quantum systems are generally formulated in terms of 
a Markovian Lindblad master equation that describes the time 
evolution of the density matrix |[T0| . 

In this article we want to investigate the interplay between 
stochastic transport processes and coherent transport present 
in the same system. Here we consider only the steady-state 
properties of the system. To treat both processes on equal 
footing, we incorporate the classical hopping terms into the 
quantum master equation. The stochastic hopping is mod- 
eled by appropriately chosen quantum jump operators. Such a 
construction has also been used to find quantum master equa- 
tions that describe a quantization of kinetic Ising models ifTTIl . 
These models obey detailed balance and allow for an exact 
solution. Considering hopping models in this more general 
quantum framework allows now for additional quantum trans- 
port, so to speak, on top of the classical hopping evolution. 
We can choose an arbitary, particle number conserving Hamil- 
tonian to mediate the coherent transport and investigate the 
effect this quantum pertubation has on the classical hopping 
process. 

Our article is organized as follows: First, in section [TTJ we 
introduce the hopping model and formulate the problem as a 
quantum master equation. We choose in different instances ei- 
ther the ASEP or the SEP as the underlying stochastic process. 



In the following section III the quantum analog of the ASEP 
is treated numerically in the framework of matrix product den- 
sity operators (MPDO). The master equation for a chain of 
N = 40 sites is evolved in time, until the steady state is 
reached. The current, the particle density, as well as the par- 
ticle density-density correlations are computed. In section IV 
the SEP is considererd. The two-point correlation functions 
of the SEP can be calculated exactly in the steady state, and 
the scaling of the current for larger lattice sizes is investigated. 
The conclusion are then drawn in the final section [V] 
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FIG. 1: spin- 1/2 chain with jump-operators and coherent evolution. 
The inflow is given by rate a, the outflow with rate f3, stochastic 
hopping to the left with rate (fL and to the right (fR. The strength of 
the coherent evolution is given by A 



II. FORMULATION OF THE QUANTUM MASTER 
EQUATION 

We study a system of hard-core particles in a one dimen- 
sional chain of length TV, where each site (1 < k < N) can 
be either occupied or empty. In the spin chain picture this cor- 
responds to either spin up (occupied), or spin down (empty). 
At the boundary k = 1, we allow for an inflow with a rate 
a and at k = N for an outflow of particles, given by a rate 
p. The particles at each site are allowed to hop stochastically 
to the left with a rate cfL and to the right with a rate (fR. In 
this article we consider only two cases. First, the fully asym- 
metric case (ASEP), where the particles only hop to the right, 
(Pr = 1 and (fL = 0. Second, the fully symmetric case (SEP) 
where both hopping rates are equal (fL = = 1- 

These stochastic particle jumps correspond to the classical 
stochastic exclusion process and are formulated in terms of 
quantum-jump operators L M . These operators govern the in- 
coherent evolution of the quantum master equation. These 
jumps can be formulated in terms of spin-flip operations: 




L[ = \fcX(7 (1) 
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Here, the a ± correspond to the Pauli raising and lowering op- 
erators. The master equation written with only these operators 
reproduces exactly the classical stochastic behavior with the 
appropriate jump rates. This can be seen, when one restricts 
oneself to density matrices diagonal in the computational ba- 
sis. With the chosen definition of the Lindblad operators, the 
classical master equation for the probabilities is reproduced. 
In this generalized framework, we can now also allow for an 
additional coherent evolution of the system, when choosing 
an appropriate Hamiltonian. The XX-Hamiltonian 
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gives rise to the free coherent evolution of the hardcore par- 
ticles. Furthermore it has the property that it conserves the 
number of particles. The full quantum master equation for the 
density operator p can be written as: 

d t p = -i [p, H} + J2 L, P L\ - \ {Lj, L„; p) (3) 

/I 

EE C[p]. 



The central observables are the particle density n k = ^X G k 
and the current j k of particles. To find the right expression 
for the particle current, we consider the continuity equation 
for the density n k . The continuity equation is obtained in the 
Heisenberg picture when the adjoined of C is acting on n k . 

d t n k = £t [n k ] (4) 

= -i[H,n k ]^^L^n k L^-]^{L^L^n k } 
\± 

Let us first consider only the (ASEP) with a single Lindblad 
operator in the bulk, since <p L = 0. We find, that n k obeys a 
continuity equation of the form: 
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We may now interpret the sum of the terms 
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as the total current-density jk,k+i of the system. Note, that 
there are two different contributions to the current, the coher- 
ent part j co due to the dynamics generated by the Hamilto- 
nian and the stochastic contribution j st originating from the 
quantum-jump operator induced hopping. We observe, that 
the stochastic contribution to the current corresponds exactly 
to the current present in the classical model ll3lfT2ll. 



III. THE ASYMETRIC EXCLUSION PROCESS 

The steady state of the ASEP, without any additional quan- 
tum evolution, i.e. A = 0, can be calculated exactly lfT2l and 
its solution can be written in terms of a matrix product density 
operator (MPDO) |[T3l . The general form of a MPDO is given 

by, 
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where d denotes the dimension of the local Hilbert space and 
Mi k j k +i are D k x D k+ i dimensional matrices, In the case of 
the ASEP, the steady state solution of the master equation is 
given by the choice, see Appendix [A] 
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FIG. 2: Classical phase diagram of the asymmetric exclusion pro- 
cess. The ASEP processes three phases, LD (low-density), HD (high- 
density) and the maximal current phase MC. 



response only depends on the absolute value of A. The simu- 
lations were done for a lattice with N = 40 sites. The matrix 
bond dimension of the MPDO was chosen as D = 60 and we 
chose a Trotter step At — 10 -4 . 

To get an understanding of how the system responds to the 
quantum perturbation in each phase respectively, we have cho- 
sen to compute the steady state at different values of a and (3, 
which correspond to points chosen to lie in different phases. 
Four different points were chosen. Corresponding to the co- 
existence line (CL), where a = (3 and (3 + a < 1, we chose 
the point a = (3 = 1/4. In the maximum current phase (MC), 
we chose a = /3 = 3/4. For the low-density (LD) and the 
high-density (HD) phase, we chose a = 1/4 /3 = 1/2 and 
a = 1/2 (3 = 1/4 respectively. 

The observable we considered first was the density distribu- 
tion (rik) = Tr {rikp} as a function of the lattice site k, FIG 

13 



which is a MPDO with a matrix dimension of D = N + 1. 
In lfl~2l the classical phase diagram with respect to a and (3 
has been found FIG. [2] In general there are three phases, the 
low-density phase LD, the high density phase HD, and the 
maximal current phase MC. The three phases exhibit different 
behavior in the current and density. The statistical current and 
the corresponding density in the large N limit are given by, 
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as has been calculated in [12]. We want to understand the sys- 
tem's response to a quantum mechanical perturbation at dis- 
tinct points in the phase diagram. 

Note, there is a special line in the classical phase diagram 
marked by a + (3 = 1. Along this line the density operator 
for the classical probabilities factorizes, and mean-field the- 
ory becomes exact, Appendix [A[ If we now add the quantum 
perturbation, we see, that the state, 
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is still the steady state of the system for arbitrary values of A. 
That is, when the system is classically uncorrelated, the quan- 
tum perturbation has no effect on the system. 
To see that in the regime, where the stochastic steady state 
is correlated, the coherent evolution alters the steady state, 
we need to calculate the steady state of the system numeri- 
cally by time-evolving the density matrix, until we reach the 
steady state. The numerical simulations of the real time evo- 
lution is performed by a modification of the TEBD algorithm 
for mixed states CUE), see Appendix [B] for details. As ini- 
tial state for the evolution we have chosen the classical steady 
state |7|. We then changed the value ofA = 0to,A = l/2 
and A = 1, and evolved the MPDO, until the steady state was 
reached, i.e. until all considered observables did not change 
any more. Negative values of A were also simulated and led 
to the same results. We conclude from this, that the systems 
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FIG. 3: (color online) Density distribution (rik) for the different 
points (a, 13), (a) HD (1/2,1/4), (b) LD (1/4,1/2) and (c) MC 
(3/4, 3/4)as well as (d) CL (3/4, 3/4). The black solid line cor- 
responds to A = 0, i.e. the classical solution. The red dashed line 
corresponds to a quantum perturbation with A = 1/2, and the blue 
dashed-dotted line to a perturbation A = 1. 



Furthermore, we calculated the values of the two-point cor- 



relation functions, of the densities rik 
of sites, 
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(10) 



The expectation values are taken with respect to the system's 
steady state. In the figures FIG. |4|5|6|7| the correlation func- 
tions are compared to the different contributions to the current 



defined in for different values of A = 0, 1/2, 1. 



FIG. 4: (color online) HD phase: The left column shows the density- 
density correlations ([ToJ for a = 1/2 and ft = 1/4 (HD) and dif- 
ferent values of coherent driving A. In the right column the current- 
density is plotted as a function of the lattice site k. The black dash- 
dotted line depicts the total current j tot . The blue solid line cor- 
responds to the stochastic contribution j st and the red dashed line 
shows the coherent contribution j co . The plots are ordered from top 
to bottom as (a) A = (b)A = 0.5 and (c)A = 1. 



FIG. 5: (color online) LD phase: In the left column we plotted the 
density-density correlations (10 1 for a = 1/4 and (3 = 1/2 (LD). 
The right column depicts the current-density as a function of the lat- 
tice site k. The black dash-dotted line amounts to the total current 
j tot . The blue solid line corresponds to the stochastic contribution 
j st and the red dashed line shows the coherent contribution j co . The 
plots are ordered from top to bottom as (a) A = (b)A = 0.5 and 
(c)A = 1. 



The first observation to be made is that the individual contri- 
butions to the total current are no longer constant throughout 
the lattice anymore, they show a dependence on the lattice 
site. The total current however, i.e. j tot = j co + j st , is still 
constant at each site of the lattice, as is required, since the 
system is in a steady state. 



a. The high- and low -density phases FIG. ?J5' These 



two phases are, just as in the classical set up, related by a 
particle-hole transformation and exchanging the numbering 
of the lattice sites from left to right. All the plots reflect this 
symmetry. One observes that the correlation functions in the 
classical regime, A = 0, are already quite short ranged and de- 
cay rapidly. The quantum perturbation in both cases leads to 
a further decay of the correlations. The total current remains 
stable with respect to the quantum perturbation and does not 
change its value notably. The individual constituents to the 
current however change their behavior. At the boundaries the 
stochastic contributions are increased, whereas the coherent 
current gives rise to a flow in the opposite direction. 



b. The maximal current phase FIG. [6j In the classical 
process, the MC phase corresponds to the maximum amount 
of current the system can carry. Allowing for a quantum per- 
turbation, the system makes use of the additional transport 
capacity and increases its total current. For these boundary 
conditions the stochastic as well as the coherent contributions 
flow in the same direction. Note, however, that the system 
in the bulk prefers to use the stochastic transport, whereas at 
the boundaries coherent transport is preferred. The classical 
correlation function initially assumes negative values close to 
the boundaries. The onset of the quantum perturbation also 
reduces the magnitude of the correlations in this phase, even 
though the total amount of current is increased. 

c. Coexistence line FIG. [7[ For the chosen boundary 
conditions, the amount of correlations initially present in the 
steady- state are decreased, when switching on the quantum 
perturbation. The final steady-state, were A = 1, however, 
still shows the presence of correlation to a higher degree than 
in the other phases. The total amount of current carried by the 
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FIG. 6: (color online) MC phase: The left column shows the density- 
density correlations ([ToJ for a = 3/4 and /3 = 3/4 (MC) and dif- 
ferent values of coherent driving A. In the right column the current- 
density is plotted as a function of the lattice site k. The black dash- 
dotted line depicts the total current j tot . The blue solid line cor- 
responds to the stochastic contribution j st and the red dashed line 
shows the coherent contribution j co . We plotted three different val- 
ues for (a) A = (b)A = 0.5 and (c)A = 1. 



system, however, is decreased. The coherent contribution is 
negative throughout the system. 



IV. THE SYMMETRIC EXCLUSION PROCESS 



If we allow for stochastic hopping in both directions with an 
equal rate cp and turn off the coherent evolution, the model de- 
scribes the classical symmetric exclusion process. The sym- 
metric exclusion process is known to possess only a single 
phase classically with a vanishing current |5 |. We now con- 
sider the quantum perturbation to the system. Note, that 
now the quantum-jump operators are related via L kk+1 — 

L L \ k+1 = Lj^fc+i. This allows us to rewrite the full master 




FIG. 7: (color online) Coexistence line: The left column shows the 
density-density correlations |To| for a — 1/4 and f3 = 1/4 (CL). In 
the right column the current-density is plotted. The black dash-dotted 
line amounts to the total current j tot . The blue solid line corresponds 
to the stochastic contribution j st and the red dashed line shows the 
coherent contribution j co . The values for the coherent driving are 
chosen as (a) A = (b)A = 0.5 and (c)A = 1. 

equation as, 

d tP = -i[ P ;H] (11) 

N-l 

+ Yl \.l L k,k+u p]\ £t,fc + i] + [[L ] Kk+1 ]p]\L Kk+1 ] 
k=l 

+ L x pL\ - ^{L{, Li; p} + L N pL ] N - ]^{L ] N , L N ; p}. 

It is possible to calculate the nearest neighbor two-point cor- 
relation functions exactly. To see why this is possible, we first 
transform the Pauli raising and lowering operators, a + and 
a to fermionic modes by means of the Jordan- Wigner trans- 
formation. The fermionic modes read then, 
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One can verify, that these modes now obey the fermionic 
anti-commutation relations, {dk',a\} = Skj. It is possible 
to calculate the evolution of the fermionic two-point function 
( a k a m) from the master equation via, 

d t {a\a m ) = Tr {tf[a\a m ]p } = (tf[a{a m ]). (13) 

Since the commutator of two pairs of fermionic modes is 
again an operator made up from two fermionic modes, we see, 
that the time-evolution of the fermionic two-point functions 
again only depends on two-point functions. So the two-point 
correlation functions of the steady state can be computed ex- 
actly. The equations for the correlation functions read, 
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$l,N)(atn a l) 



{[(a, 



m+ ia m+1 



+1 a m+ i/ 



2 

+ ( a m-l^m-l) _ ( a ln-l a rn-l}}3l,r, 
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(14) 
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+iA ((a^a^) + (a^aj^) + (<4_ 1 a z t ) + (4i+i<4>) 15 ) 

The other correlation functions are related to these two by the 
identities imposed due to the anti-commutation relations of 
the fermionic modes, thus (a^aj) = (a m ai)* and (a^ai) = 

The steady- state correlations can be computed from these 
equations by requiring stationarity, i.e. d^a^ai) = 
dtia^aj) = 0. One needs to solve the corresponding sys- 
tem of linear difference equations. 

The current as well as the particle number density can be ex- 
pressed in terms of these correlators. One finds for the particle 
number density (n^) = (a\ak) and the two contributions to 
the current read, 



ASt 

Jk,k+1 



(16) 



= cp((n k ) - (nfe+i)) 

3k%+i = x (<4 a *+i> - w4+i)) • 

Note, that the stochastic current now only depends on the dif- 



ference of the densities at adjacent sites. The equations (14) 



and ( p3) are solved numerically for different parameters a and 

The classical SEP without any further driving obeys the de- 
tailed balance condition, and thus does not support a steady 



state current. When one allows for an external driving of the 
particles at the boundaries, as we do in our example, a current 
is induced in the SEP steady state. This current however van- 
ishes as ~ 1/N in the system size N, see FIG.[8ja) . If we turn 
to the coherent part, we also see this behavior, i.e. the current 
vanishes in the limit N — >> oo FIG.[5Jb). We deduce from this, 
that the quantum perturbation to the SEP is an irrelevant per- 
turbation and does not lead to a qualitatively different behav- 
ior of the systems transport properties. This point is also af- 
firmed by considering the particle density distribution as seen 
in FIG.|9fa). Here one sees that the density distribution in the 
presence of a quantum perturbation almost coincides with the 
classical distribution. The different contributions to the cur- 
rent, FIG.|9jb), are identical in the system's bulk, and add up 
to an increased current J tot . Whether the quantum perturba- 
tion is completely irrelevant however can not be deduced from 
just considering the steady state density and the current alone. 
Here one would need to take higher order correlations into ac- 
count, as for instance the current-current correlation function 
at unequal times. 



(b); 



FIG. 8: (color online) The plots demonstrate the scaling of the cur- 
rents for large N. Plot (a) depicts the scaling of the total current 
jtot _ jst ^uQ dashed line for A = 0. The plot (b) corresponds to 
A = 1. Here the total current j tot is represented by the black dash- 
dotted line. The coherent contribution corresponds to the blue solid 
line. The boundary values were chosen asa = /3 = 3/4 



V. CONCLUSIONS AND OUTLOOK 

We have investigated a quantum perturbation to the dynam- 
ics of the stochastic asymetric exclusion process as well as 
to the symmetric exclusion process. We find that we can 
rephrase the stochastic master equation as a quantum equa- 
tion that fully reproduces the classical dynamics. The quan- 
tum perturbations modify the steady state behavior and allow 
for two different types of currents, which each on their own 
can vary as a function of the site. Numerical simulations of 
the full master equation indicate, that the underlying classi- 
cal phase-diagram of the stochastic process is respected. The 
steady state responds to driving due to the boundary terms 
with a different behavior in current and density. A further step 
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FIG. 9: (color online) Plot (a) shows the density distributions for the 
values A = 0, 1 and for a = f3 = 3/4 . Calculations for larger sys- 
tem sizes indicate, that this deviation vanishes in the thermodynamic 
limit. The second plot (b) shows the dependence of the current dis- 
tribution on the lattice site for A = 1. The stochastic contribution j st 
as well as the coherent contribution j co only marginally differ at the 
boundaries 



would be to investigate the current-current correlation func- 
tion of the SEP, to see whether the quantum perturbation has 
an effect to the current-fluctuations. Furthermore, other, more 
complex models with an interplay between stochastic and co- 
herent dynamics can be investigated along these lines. 
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Appendix A: Derivation of the steady state 

The steady state can be calculated in matrix product form by 
applying the Derrida method of the quadratic algebra lfT2l . 
To better understand this, we shall rewrite the MPDO with 
translationally invariant matrices Mij as, 



(Al) 



This MPDO now serves as steady state ansatz for the master 
equation. To derive the quadratic algebra, consider the total 
master equation ([3]), where we now regrouped the terms in 
local two-body interactions and the single particle boundary 
terms. This can be written as, 



N-l 



C[p] = CM + £k,k+i\p] + £i[p]- (A2) 



fc=l 



We now try to find relations for the matrices Mij so that p 
will correspond to the steady state of the system. We introduce 



additional ancilla matrices and require, that 
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Furthermore we require, that the single site operators at the 
boundaries have to satisfy, 



(V\d 



M , M ,i 
M lj0 Mi i 



= (V\ 



'M 0>0 M ,i 
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W) = - 
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We see that the total sum ( |A2[ ) telescopes to zero and p is 
the steady state solution of the equation. To find an explicit 
solution one needs to construct a representation of the matri- 
ces that obey the given algebraic constraints. In the sim- 
ple case, where A = and cp = 1, we reproduce the known 
classical algebra for the ASEP steady state lfT2l . Here one 
chooses M ,i = Mi >0 0, as well as M ,i M li0 0. 
The diagonal terms, upon stetting, Mo,o = — Mi } i = 1, 
have to satisfy the algebra Mi 5 iMo,o = Mo,o + M^i. The 
boundary terms have to fulfill, (V | M ,o = l/a(V\ and 
Mi } i | W) = 1/^| W). One possible representation for this 
algebra is given by ([5]), when considering TV — ^ oo. Note, 
that when one wants to reproduce the steady state of a finite 
system of size N, it suffices to choose a matrix dimension of 
D = N + 1 for this specific representation. One can also ver- 
ify, that at a + f3 = 1 the representation of the algebra can be 
chosen one-dimensional |[T2ll . 



Appendix B: Numerical Implementation 

For the numerical evolution of the density-matrix according 
to ([3]), we apply a numerical scheme developed in fT3l . Start- 
ing from the initial density-matrix po given as an MPDO, we 
apply the CP-map £(t) = exp(tC) for a small time step At 
and approximate the resulting density operator, that has now 
an increased bond dimension, with an MPDO that has a bond- 
dimension Dk corresponding to that of the original MPDO. 
The approximation of the operator p(t + At) = 8 {At) p(t) 
is chosen, such that the Hilbert- Schmidt norm \\p(t + At) — 



-'new \\HS 



Tr 



is minimized.This 



(p(* + A*) 

optimization can be performed efficiently by sweeping from 
left to right over the individual sites and optimizing the ma- 
trices locally. For the application of the CP-map to be com- 
putable, we perform a second-order Trotter expansion of the 
CP-map as follows, 

£(£, At) ~ £(C , At/2)£(£ e , At)£(C Q , At/2), (Bl) 

where C = C e + C Q corresponds to splitting the Liouvillian 
into commuting terms which act on the sites (2fc,2fc + l) and 
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(2k — 1, 2k), respectively. The resulting MPDO p new is then is repeated, 
chosen as initial condition for the next step and the procedure 
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